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We predict strong non-local effects in the three-terminal hybrid device, comprising the quantum 
dot embedded between two conducting leads and third superconducting reservoir. They result from 
competition between the ballistic electron transfer and the crossed Andreev scattering. The non-local 
voltage induced in response to the ’driving’ current changes the magnitude and sign upon varying 
the gate potential and/or coupling to the superconducting lead. Such effect is robust both in the 
linear and non-linear regimes, where the screening and the long-range interactions play significant 
role. This novel subgap transport is provided by the Shiba states and can be contrasted with much 
weaker non-local effects observed hitherto in the three-terminal ’planar’ junctions. 


Multi-terminal systems enable measurements of both the local and the nonlocal voltages/currents between selected 
electrode pairs [lj . The non-local transport of charge 00, heat 0 and spin jff] via hybrid devices consisting of the 
normal and superconducting reservoirs are currently of interest for the basic research and innovative applications. 
Electrons traversing metal-superconductor interface are glued into the Cooper pairs, and conversely, the Cooper pairs 
are split into the individual electrons 0- In both processes there emerge the entangled carriers, leading to nonlocal 
correlations. These effects can be amplified by inserting the quantum dots between the reservoirs [ill] . In this reg ard, 
the three-terminal structures are especially useful, because they allow for efficient splitting of the Cooper pairs [12|- 14]. 
give rise to spin filtering [15,j , generate the correlated spin currents fill. se parate the charge from heat currents 17j , 
enable realization of the exotic Weyl or Majorana-type quasi-particles jig], etc. 

Very spectacular non-local effects are provided by the crossed Andreev reflections (CAR), operating in a subgap 
regime. The ’driving’ current applied to one side of the multi-terminal junction can yield either positive or negative 
nonlocal voltage response at the other interface, depending on a competition between the ballistic electron transfer 
(ET) and the CAR processes. Such changeover has been observed in three-terminal planar junctions , using a 
piece of superconducting sample sandwiched between two conducting (normal or magnetic) electrodes. The induced 
non-local conductance, however, was much weaker from the local one in agreement with theoretical predictions [l9i - l2li ]. 



FIG. 1: (color online) Scheme of the three-terminal device consisting of two conducting leads (L and R), superconducting 
reservoir ( S ) and the quantum dot ( QD ). The ’driving’ current in the L — QD — S loop induces the non-local voltage ’response’ 
of the floating R electrode. 

Here we propose a different configuration, where the quantum dot is built into the three-terminal hybrid as sketched 
in Fig. [I] Proximity effect converts the quantum dot into,a kind of, superconducting grain and its subgap spectrum 
develops the, so called, Andreev or Shiba bound states [22], which substantially enhance the non-local transport. 
We show that effective non-local conductance can be comparable to the local one and can change sign from the 
positive to negative values by increasing the coupling Tg to superconducting electrode or by appropriate tuning of 
the gate potential. The gate potential is also controlling symmetry of this effect. Experimental tests of such effects 
should be feasible using the three-terminal architecture with such quantum dots as the carbon nano-tubes ]7l. [23j, 
semiconducting nano-wires 0Ei] or self-assembled InAs islands ]26]. 
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Microscopic model 

Some aspects of the local and non-local transport properties for this three-terminal device could be inferred by 
extending the Landauer-Biittiker approach [27H3lj (see the 1-st subsection of Methods). On a microscopic level, we 
describe this system in the tunneling approximation [32] by the Hamiltonian 

H = 'y ] £qkc]xktT C <xkcr + ^ [ [ e 0 ~ HJ {v)\ <$ a d a + ^ [ (^ac] yka d a + t a dl T C a ka-'^ 

<y—L,R k,cr cr a,k,a 

+ ^SkC^ ska Cska - ^2 ^ ( C S-fcT C S'fc! + C SklCS-k2j (1) 

k,cr k 

with standard notation for the annihilation (creation) operators of the itinerant c^2 and localized dot d electrons. 
The first term describes the left ( a = L) and the right ( a = R) conducting leads. The subsequent term refers to 
the quantum dot (QD ) with its energy level eo shifted by the long-range potential U( r). Hybridization between the 
QD and itinerant electrons is characterized by the matrix elements t a . The last two terms in m correspond to 
the BCS-type superconducting reservoir with an isotropic energy gap A. Addressing here the subgap (low-energy) 
transport we assume the constant tunneling rates T a = 27r \t a \ 2 5(E — e ak ) = 27r|t a | 2 p a , where p a is the (normal 
state) density of states of a lead. In what follows, we assume the superconducting gap A to be the largest energy 
scale in the problem. 


Subgap charge transport 


The charge current J Q flowing from an arbitrary lead a = { L , R , S} can be evaluated using the Heisenberg equation 
J a — e(N a ) = —ie/h{[N a , H]) [33]. In particular, the current Jl(r) from the normal L ( R ) electrode is given by .33] 


J HR) = f 


■HR)' 


f s~ir 
JL(R)^ 1 


I G< 
2 Gl1 


( 2 ) 


where and denote the matrix elements (in the Nambu representation) of the retarded and lesser QD Green 
functions, respectively. From now onwards we consider the current Jl focusing on the subgap voltage, smaller than 
the energy gap |A|. In such regime there are possible: the ballistic electron transfer (ET) from L to R electrode, the 
direct Andreev reflection (DAR) when electron from L lead is converted into the Cooper pair in S reservoir and hole 
is scattered to L electrode, and the crossed Andreev reflection (CAR) which is similar to DAR except that hole is 


scattered to R electrode. They can be expressed as [34 [ 


J E T = f / dET L T R \ G ll \ 2 (fL - Ir) , 

J D A R = 2 e_J dErllGTnl^h-h), 

jC A r = %J dE y l T r \ G \ 2 \ 2 (/i-/ fl ) , 


( 3 ) 

( 4 ) 

( 5 ) 


where f a = f a (E) = {exp[(£- n a )/k B T] + 1} 1 and f a = f a (E) = 1 - f a (-E) = {exp[(E + p a )/k B T} + 1} 1 are the 
Fermi-Dirac distribution functions for electrons and holes, respectively. Let us remark, that only the ET |3]) and CAR 
(fo|) contributions lead to the non-local effects, because they depend on the chemical potentials of both conducting 
(L and R) electrodes. Since these ET and CAR processes deliver different types of the charge carriers to the right 
electrode, the induced voltage V B would be a probe of the dominant transport mechanism. 

Relationship between the ET and CAR processes can be inspected by studying their transmissions, defined as 
T et (E ) = Fir^|Gi 1 (E)| 2 and T CAR (E) = TlTr\G\ 2 {E)\ 2 (see Fig. O. Deep in a subgap regime (i.e. for \E\ <C |A|) 
the Green function G r (E) simplifies to the familiar BCS structure [13 ■ Its diagonal part is given by G r u (E) = 
u 2 / [E - Ea + §r,/v] + v 2 / [E + Ea + §rJv] with the quasi-particle energy Ea = + A q D , where A qd = ^Tg. 

Subgap spectrum consists thus of two Shiba states at ±Ea whose spectral weights are u 2 = \ [1 + A qd/Ea] and 
v 2 = 1 —u 2 with the quasiparticle broadening T n = El+F/j. The single electron transmission T et {E) is a quantitative 
measure of this subgap spectrum. The left panel in Fig. [2] illustrates evolution of the Shiba states upon increasing 
the coupling Tg while the right panel shows a transfer of the corresponding spectral weights u 2 -O- v 2 upon changing 
the QD level eo by an applied gate voltage. 
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FIG. 2: (color online) Transmissions of the ET (dashed lines) and CAR (solid lines) transport channels obtained at zero 
temperature for Tr = Tl- The left panel refers to eo = 0 and the right one to Fs = 3 Fl. 


Transmission of the anomalous CAR channel, on the other hand, depends on the off-diagonal part of the matrix 
Green function G\ 2 (E) = uv/ [E — Ea + fTjv] — uv/ [E + Ea + |Tjv] , where uv = \ Aqd/Ea■ It also has maxima 
around the same Shiba states EE a but with a different amplitude, sensitive to the induced pairing (dj,d-[-). This is 
a reason why T CAR (E) quickly diminishes whenever Tg is decreased or the QD level eo departs from /xg = 0 (solid 
lines in Fig. H. 

Confronting both these transmissions reveals that the non-local transport predominantly comes from the CAR 
process when the coupling Tg (to superconducting electrode) is sufficiently strong and the QD level eo is close to 
the chemical potential /xg. Otherwise, the non-local effects are dominated by the single electron tunneling (ET). The 
related changeover can be detected by measuring the voltage Vr in the floating R electrode, in response to the current 
in the L — QD — S branch. Such voltage Vr can vary between the positive and negative values and the non-local 
resistance can be tuned by the gate potential lifting/lowering the Shiba energies. 


Linear response 

Practical realizations of the setup (Fig. [T| would allow to measure the local and the non-local resis¬ 
tances/conductances within the four-probe scheme I27| - f30| . where the potentials and currents are treated on equal 
footing (see the Method). In a weak perturbation limit the response would be linear 

Jl * eC E L T R {V L - V R ) + eC EAR [{V L - V s ) - {V s - V L )\ + eC c L AR [{V L - V s ) - (Eg - Vr)] . (6) 

The coefficients for /3 = ET, DAR or CAR can be determined from the equations (131151) and they read 




s 0 /r L 


FIG. 3: (color online) The non-local resistance 2 Rrs,ls / Dr as a function of Ts/Fn ratio (left panel) and the QD dot level eo 
(right panel) obtained in the linear limit for three representative temperatures. 
























4 



(7) 


At zero temperature —^ ~ 5(F), hence depend on the transmissions T^{E —> 0). 

Treating the potential Vg as a reference level we analyze the induced voltage Vr in response to the ’driving’ 
current Jl = Jls- The local resistance (Vl — Vs)/ Jls = Rls,ls is due to the DAR processes whereas the non¬ 
local one [Vr — Vs)/ Jls = Rrs,ls results from the single electron tunneling (ET) competing with the anomalous 
crossed Andreev reflection (CAR) processes. Fig. [3] shows this non-local resistance Rrs,ls normalized with respect 
to D r /2 = Rls,lrRrs,ls + Rls,lrRrs,rl + Rrs,lsRrs,rl [defined by equation (pj) in Methods]. The left panel 
shows that Rrs,ls has a negative sign (signifying the dominant CAR processes) only for sufficiently strong coupling 
Tg > Tat. This is a straightforward consequence of the (zero-energy) ET and CAR transmissions (Fig. [2]). The right 
panel of Fig. [3] displays the non-local resistance versus the QD level eo- In the linear regime the negative nonlocal 
resistance occurs when eo ~ Ms for sufficiently strong coupling Tg > Tat. Since Tg and eo can be experimentally 
varied in the realizations of the superconducting-metallic devices with the quantum dots [Tj, [23|i26| , such qualitative 
changes should be observable. 


Beyond the linear response limit 

To confront these findings with the non-local effects observed so far in the ’planar’ junctions we now go 

beyond the linear response framework. For arbitrary value of the ’driving’ voltage Vl we computed sclf-consistently 
Vr , guaranteeing the net current Jr to vanish. Under such non-equilibrium conditions the long-range potential U( r) 
plays an important role in the transport when the charges pile up in the electrodes and the quantum do t [368 . It 
affects the chemical potentials and the injectivities of the leads and contributes to the screening effect [37l - l4(j ]. The 
potential U(r) has to be properly adjusted, depending on specific polarization of the system [38j] (for details see the 
2-nd subsection of Methods). 




eV L / r L eV L / r L 

FIG. 4: (color online) The non-local voltage Vr (left panel) and its derivative dVa/dVL (right panel) induced in the floating R 
lead in response to the ’driving’ voltage Vl- 

Figure [4] shows the induced non-local voltage Vr and its derivative with respect to Vl for several couplings Tg 
and temperatures, obtained for U( r) = 0. At low voltage \Vl\ the induced potential Vr is proportional to Vl, as 
we discussed in the linear response regime (Fig. [3]). Upon increasing the ’driving’ voltage |Vj,| the Shiba states EE a 
(indicated by vertical lines in Fig. are gradually activated, amplifying the non-local processes. For Tg > Tat 
we hence observe local minima (maxima) of Vr at the quasiparticle energies Ea [—Ea)- Further increase of \Vl\ 
leads to revival of the dominant ET channel. The derivative dVg/dVL, which is related to the ratio of the local and 
non-local differential resistances Rls,rs/Rls,ls, can be measured by the standard lock-in method. Our results differ 
qualitatively from the properties of the planar junctions (where the ET and CAR dominated regions are completely 
interchanged) UHl because the non-local transport occurs through the Andreev states, that are localized at two 
normal-superconductor interfaces separated by a distance d comparable to the coherence length of superconductor. In 
consequence, the anomalous CAR transport is possible only for eVg exceeding the characteristic Thouless energy 0- 
0 . 
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eV L / r L eV L / r L 

FIG . 5 : (color online) The non-local voltage Vr and its derivative with respect to Vl obtained at low temperature for eo = 0 
(left panel) taking into account the screening effects U{ r ). The lower panel shows dVR/dVL for eo/F l = ± 1 . 


Feedback effect of the long-range potential U(r) = U eq + u a F a (where U eq denotes the equilibrium value 
incorporated into eo) is illustrated in Fig. [5] The quantitative changes are observed for all voltages, however, the 
qualitative behaviour is similar to that found in the linear regime (Fig. [4]). The screening effects and injectivities are 
calculated here in the self-consistent way (di^-HTif (discussed in the 2-nd subsection of Methods). This selfconsistent 
treatment of U( r) partly suppresses both the non-local voltage Vr and dVR/dVL■ The right panel of Fig. [5] shows 
dVR/dVL with respect to Vl outside the particle-hole symmetry point, i.e. for eo = ±Ti. These asymmetric curves 
can be practically obtained by applying the gate potential to the quantum dot. 


Summary and Outlook 


We proposed the three-terminal hybrid device, where the quantum dot is tunnel-coupled to two normal and another 
superconducting electrode, for implementation of the efficient non-local transport properties. We investigated such 
effects in the linear and non-linear regimes. We found that in the both cases the non-local resistance/conductance 
can change from the positive (dominated by the usual electron transfer) to negative values (dominated by the crossed 
Andreev reflections) upon varying the coupling to superconducting electrode Tg and tuning the QD level eo- 

This nano-device would enable realization of the strong non-local conductance (comparable to the local one) by 
activating the Shiba states formed at sub-gap energies ±Ea- They substantially enhance all the transport chan¬ 
nels, in particular promoting the CAR mechanism (manifested by the negative non-local conductance/resistance) 
when the coupling to superconducting electrode is strong Tg > Tg + T#. We predict the negative non-local con¬ 
ductance/resistance both, in the linear regime and beyond it. For the latter case such behavior would be observable 
exclusively in the low bias voltage regime \Vl\ < Ea/& capturing the Shiba states. The quantum dot level eo (tunable 
by the gate potential) can additionally control asymmetry of the non-linear transport properties, affecting the CAR 

transmission T gar (±Ea ) oc 1 + (2eo/r $Y 


Strong non-local properties of the nano-device (shown in figure 1) can be contrasted with the previous experimental 
measurements for the three-terminal planar junctions (consisting of two N — S interfaces separated by a supercon¬ 
ducting mesoscopic island) @41. Russo et al. Q reported evolution from the positive to negative non-local voltage Vr 
induced in response to the ’driving’ bias Vl- At low Vl the ET processes dominated, whereas for higher Vl the CAR 
took over. The sign change of Vr occurred at voltage Vl related to the Thouless energy (such changeover completely 
disappeared when a width of the tunneling region via the superconducting sample exceeded the coherence length). 
Similar weak negative non-local resistance/conductance has been observed in the spin valve configurations [I, [a]. In 
the planar junctions the non-local conductance was roughly 2 orders of the magnitude weaker than the local one Q. 

Summarizing, we proposed the nanoscopic three-terminal device for the tunable (controllable) and very efficient 
non-local conductance/resistance ranging between the positive to negative values. Our theoretical predictions can 
be verified experimentally (in the linear response regime and beyond it) using any quantum dots Fri. [23142(1 attached 
between one superconducting and two metallic reservoirs. Such measurements are called for. 

























6 


Appendix A: Landauer-Biittiker formalism 


The four-point method 27] is well established technique for measuring the resistance in a ballistic regime. Voltage 
Vki measured between k and l electrodes in response to the current Jy between i and j electrodes defines the local 
(ij = kl) or non-local (ij ^ kl ) resistance via 


Ri.j.kl — 


_ Vk — Vi fik ~ Hi _ ^Hki 


Ji 


e J,;- 


cJij 


( 8 ) 


where A/ifcj = /jj - /j; is a difference between the chemical potentials of k and l electrodes. The formalism has 
been later extended by Lambert et al. p8l . [20t to systems, where electron tunneling occurs between one or more 
superconductors. The current from i-th lead depends on the chemical potential hs of superconducting reservoir, 
because the scattering region acts as a source or sink of quasi-particle charge due to the Andreev reflection (see e.g. 
Ref. [H). 

Adopting this approach, we analyze here the local and non-local transport properties of the three-terminal hybrid 
system consisting of two normal (L and R) leads coupled through the quantum dot with another superconducting (S) 
electrode. We consider the charge transport driven by small (subgap) voltages eVu = A Hkl = Hk — Hi "C A, when the 
single electron transfer to the superconductor is prohibited. In this limit the net current flowing from the normal L 
electrode consists of the following three contributions 

Jl = C-lfl (Vl - Hr) + £ll R [(hl - Hs) ~ (hs - Hl)] + £lr R [{hl ~ Hs) - {hs - Hr)] ■ ( 9 ) 

The linear coefficient Clr refers to the processes transferring single electrons between metallic L and R leads. We call 
this process as the electron transfer (ET). The other term with corresponds to the direct Andreev reflection, 

when electron from the normal L lead is converted into the Cooper pair (in S electrode) reflecting a hole back to 
the same lead L. The last coefficient describes the non-local crossed Andreev reflection, involving all three 

electrodes when a hole is reflected to the second R lead. In the subgap regime the competing ET and CAR channels 
are responsible for the non-local transport properties. 

In the same way as © one can express the current Jr. By symmetry reasons we have = C*[r and £rr R = 
£rr R , whereas the charge conservation (Kirchoff’s law) implies Js = —Jl~Jr- From these linear response expressions 
one can estimate the relevant local and non-local resistances ©, assuming arbitrary configurations of the applied 
currents and induced voltages. Experimental measurements of such resistances © can be done, treating one of the 
electrodes as a voltage probe. In our three-terminal device with the quantum dot we can assume either the metallic 
or superconducting electrode to be floating. We now briefly discuss both such options. 


Floating metallic electrode 


We assume that the superconducting lead S is grounded and treat the metallic electrode (say L) as a voltage probe. 
This means that the net current vanishes Jl = 0 and, from the charge conservation, one finds Jr = —Js = Jrs- I n 
the linear response regime © implies the following potential differences 

( 10 ) 

( 11 ) 

( 12 ) 

with a common denominator 




zJrs 
A Hls 

cJrs 

^Hrs 

eJRS 


= Rrs,rl = 


— Rrs.rs = 


rDAR , rCAR 
*~LL +L r.B 


LR 


eD 


= Rrs,ls — 


rET _ rCAR 


2 eD 


rET 
^LR 


2 C%£ R + C 


CAR 

LR 


2 eD 


A Hls A-Hrl 
eJRs eJRS 


= Rrs.rl + Rrs,ls 


D — C ET (C 

u — l -LR\ L ' 


LL 


2 c c l £ r 


r DAR\ 
^RR ) 


rCAR 

'-'LR 


(£ll r 


+ £ 


DAR\ 
RR ) 


2 Cll R 


£ 


DAR 

RR 


(13) 


According to the definition © and using (11011121) we obtain the local ( Rrs,rs) and non-local ( Rrs,rl , Rrs,ls) 
resistances for the floating L lead. Let us notice, that a sign of the non-local resistance Rrs,ls depends on a 
competition between the normal electron transfer (ET) and the crossed Andreev reflections (CAR). The local resistance 
Rrs,rs is in turn a sum of the non-local resistances Rrs.rl and Rrs,ls- For the configuration, where the other ( R ) 
metallic lead is floating we obtain the equations similar to (11011121) with the exchanged indices L -(->• R. 
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Floating superconducting electrode 


We encounter a bit different situation, assuming the superconducting S electrode to be floating (i.e. Js = 0). The 
charge conservation Jr = — Jr = Jrr and Eq. © imply 


Afj, s 


ls 


eJLR 




SR - 


s-Jlr 


tDAR , r CAR 

n n RR ~r E LR n 

-Cl LR,LS — -p;- — li-LS,LR , 

eD 

tDAR , tCAR 

d E LL "r E LR D 

rt-LR,SR ~ -—- — rtRS,RL , 


eD 




LR 


cJlr 


= Rlr.lr = 


L°£ r + 2L c l £ r + L°£ r Ah 


J LR 

eD 


’‘LS 


eJ, 


LR 


^■Fsr 

eJLR 


(14) 

(15) 

= Rlr,ls + Rlr,sr = Rls,lr + Rrs,rl -(16) 


We notice some analogy between the resistances (11411161) and the previous expressions (11011121) . The significant difference 
appears between the non-local resistances Rrs,ls mu and Rlr,sr m Because of a minus sign in (fill) the former 
configuration seems to be more sensitive for probing the local versus non-local transport properties. 


Remarks on the determination of partial conductances 


Measurements of the local/non-local resistances provide information about the competition between various tun¬ 
neling processes. Similar information can be also deduced about the linear coefficients d'l :j . Let’s combine the results 
obtained for L (or R) and S floating electrodes. We have three independent equations, but we have to determine four 
coefficients 

nDAR | nET _ 

^LL ^LR ~ 


rDAR | rET _ 

’-'RR ’-'LR ~ 


r CAR r ET _ 

^LR ~ ’-'LR ~ 

In general, we thus cannot obtain a complete information about all conductances from the separate measurements of 
the currents and voltages. This situation differs from the case when the quantum dot is coupled to all three normal 
electrodes, where electrical transport can be characterized only by three conductances. 

Fortunately, for the case with asymmetric couplings Tr r the measurements can unambiguously determine the 
partial conductances 


Rlr,rs — 2Rrs,ls 

Ld~ r ’ 

Rlr,ls + 2 Rrs,ls 


eD r 


2 Rrs,ls 
eD R 


(17) 
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ET 

LR 


= P r ET 

— e>- L R 


Q 


DAR 

LL 


= eL 


DAR 

LL 
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DAR 

RR 


= eC 


DAR 

RR 


G 


CAR 

LR 


= eC 


CAR 

LR 


r 2 l {Rlr,ls + 2Rrs,ls) + r 2 r {Rlr,rs - 2Rrs,ls ) 


r 2 l (Rlr, 


LS 


(rl-r 2 R )DR 

Rlr,rs) 


(T| -T* r )Dr ’ 

r 2 

1 R r>DAR 
Fa'dLL > 

1 L 

r 2 l Rlr,ls + r rRlr,rs 
01-r 2 r )Dr 


(18) 

(19) 

( 20 ) 
( 21 ) 


Some inconvenience is related to the fact the tunneling rates Tx,, T# must be measured as well. 


Appendix B: Non-linear transport 

The non-linear effects are of vital importance in the transport studies of nanostructures inter alia due to limited 
screening of charge and access to far from equilibrium states of the system. Non-equilibrium transport driven by the 
voltage Vl (beyond the linear regime) in nanostructures is accompanied by substantial redistribution of the charges. 
This affects the occupancy of the quantum dot and leads to piling up of the charge in the electrodes. By long range 
Coulomb interactions the charge redistributions backreact on the transport properties. We shall address this effect in 
some detail. 




















Let’s note that we are considering here the charge transport driven by voltages safely below the superconducting 
gap e\V\ < |A| (practically we assume |A| ~ 10014)- Nevertheless, even at such small voltage (of the order of a 
few r L ) the pile-up of electric charges in the electrodes and the dot affects the transport by shifting the chemical 
potentials and screening the charge on the dot. This is taken into account in the Hamiltonian (jT]) by the term eU{ r). 

The effect has been considered first in mesoscopic normal systems by Altshuler and Khmelnitskii |36| , Biittiker with 
coworkers [37, IS and others !39j. It has been also explored in the metal-superconductor (two-terminal) junctions flo j. 
Here we follow |40j, assuming that the long range interactions modify the on-dot energy eo changing it to eo — et/(r). 
In equilibrium the potential U( r) has a constant value, which we denote by U eq . In the presence of the applied voltages 
14 (where a = L,R, S ) the deviations SU = f/(r) — U eq , in the lowest order, would be a linear function 



( 22 ) 


where (...)o denotes the derivative with all voltages set to zero and the gauge invariance implies that (Jf~) = 

1 1-32] • Our treatment here relies on the mean field like approximation. In the three terminal device with the quantum 
dot the single electron transport occurs between the left and right normal electrodes, while the (direct and crossed) 
Andreev processes involve the normal and superconducting electrodes. The currents d3j SI El) and the quantum dot 

charge n = 2 / || \G r n \ 1 2 (T L f L + T R f R ) + \G r 12 \ 2 {T L f L + ^r/r) depend on the screening potential C/(r). During 


the flow of carriers the deviations of 5U from the equilibrium value U eq can be related to the change of the charge 
carriers Sn by the capacitance equation Sn = CSU, where C is capacity of the system. The charge density as well 
as all currents depend on the voltages and 5U. This allows to write the relation between Sn = n — n eq , where n eq 
denotes the equilibrium (i.e. calculated for all voltages set to zero) value of the charge 



14 - n SU , 


(23) 


where n denotes the Lindhard function. Combining these equations we solve for known in the literature as 

the characteristic potentials and conveniently denoted by u a . They describe the response of the system to the applied 
voltages. One finds 


1 ( dn \ 

Ua ~ c + n \dV a ) 0 


(24) 


For the analysis of voltages induced in the R electrode as a result of current flowing in the L — S branch of the 
system we need both ul and ur. As in the earlier work jioj we assume C = 0 in the following. The inspection of 
the formula for n reveals that for the symmetric coupling T £ = T# the functions of both electrodes take on the same 
value ul = u r. The characteristic potentials enter the expression for the Green functions and as a result modify the 
relation shown in the figure [J] The modification is especially severe for Vl > T^. 

Let us note that n = is obtained from matrix elements G(-j and G\ 2 of the the Green functions as they 

depend on the potential U. The calculation of the characteristic potentials Ur/r require the derivatives of n with 
respect to voltages Vr/r, which enter the distribution functions. The characteristic functions define in turn the 
potential U = ulVl + urVr , which has to be introduced into the Green functions entering the expressions SUED 
for the currents. 
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